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Production of the P(4260) state in B meson decay 
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We calculate the branching ratio for the production of the meson Y(4260) in the decay B~ —¥ 
Y(4260) K~. We use QCD sum rules approach and we consider the Y(4260) to be a mixture 
between charmonium and exotic tetraquark, [cg][gc], states with J PC = 1 . Using the value of 

the mixing angle determined previously as: 9 = (53.0 ± 0.5)°, we get the branching ratio B(B —» 
Y(4260)A) = (1.34±0.47) x 10~ 6 , which allows us to estimate an interval on the branching fraction 
3.0 x 10 -8 < B y < 1.8 x 10” 6 in agreement with the experimental upper limit reported by Babar 
Collaboration. 


The F(4260) state was first observed by BaBar col¬ 
laboration in the e + e _ annihilation through initial state 
radiation [lj, and it was confirmed by CLEO and Belle 
collaborations [|j. The F(4260) was also observed in 
the B~ -A Y(4260)R _ -A K~ decay Q, and 

CLEO reported two additional decay channels: J/4/7r°7r° 
and J/4/AT+/W Q. The F(4260) is one of the many 
charmonium-like state, called X , Y and Z states, re¬ 
cently observed in e + e _ collisions by BaBar and Belle 
collaborations that do not fit the quarkonia interpreta¬ 
tion. The production mechanism, masses, decay widths, 
spin-parity assignments and decay modes of these states 
have been discussed in some reviews [H3- The F(4260) 
is particularly interesting because some new states have 
been identified in the decay channels of the F(4260), like 
the F+(3900). The Z+(3900) was first observed by the 
BESIII collaboration in the mass spectrum of 

the F(4260) -A decay channel Q- This struc¬ 

ture, was also observed at the same time by the Belle 
collaboration j9| and was confirmed by the authors of 
Ref. [13 using CLEO-c data. 

The decay modes of the F(4260) into J/^> and other 
charmonium states indicate the existence of a cc in 
its content. However, the attempts to classify this 
state in the charmonium spectrum have failed since the 
11/(35), 4/(2 D) and 4/(45) cc states have been assigned 
to the well established 4/(4040), 4/(4160), and 4/(4415) 
mesons respectively, and the prediction from quark mod¬ 
els for the 4/(3 D) state is 4.52 GeV. Therefore, the mass 
of the F(4260) is not consistent with any of the 1~~ cc 
states [1, @|. 

Some theoretical interpretations for the Y (4260) are: 
tetraquark state [0|, hadronic D\D, DqD* molecule 
[l2|, Xciw molecule~jl3t . XaP molecule fl4fl. J/i/;/o(980) 
molecule [ljj_, a hybrid charmonium [l6|, a charm- 
barvonium [l7j| . a cusp fl8l - l20l| . etc. Within the available 
experimental information, none of these suggestions can 
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be completely ruled out. However, there are some calcu¬ 
lations, within the QCD sum rules (QCDSR) approach 
[sL [2l| - l23[ . that can not explain the mass of the F(4260) 
supposing it to be a tetraquark state [13, or a D\D, 
D 0 D* hadronic molecule [24[, or a •7/'0/o(98O) molecular 
state [25|. 

In the framework of the QCDSR the mass and the 
decay width, in the channel J/ijiTCK, of the F(4260) were 
computed with good agreement with data, considering it 
as a mixing between two and four-quark states [13 • The 
mixing is done at the level of the hadronic currents and, 
physically, this corresponds to a fluctuation of the cc state 
where a gluon is emitted and subsequently splits into a 
light quark-antiquark pair, which lives for some time and 
behaves like a tetraquark-like state. The same approach 
was applied to the X(3872) state and good agreement 
with the data were obtained for its mass and the decay 
width into J/i/jTnr [27j, its radiative decay [13. and also 
in the X(3872) production rate in B decay [291 ] . 

In this work we will focus on the production of the 
Y(4260), using the mixed two-quark and four-quark pre¬ 
scription of Ref. [13 to perform a QCDSR analysis of 
the process B~ -a F(4260)IW. The experimental upper 
limit on the branching fraction for such a production in B 
meson decay has been reported by BaBar Collaboration 
Q, with 95% C.L., 

B Y <2.9 x 10~ 5 (1) 

where B Y = B{B~^ A'-Y(4260),F(4260)-A J/^tt+tt"). 



FIG. 1. The process for production of the Y(4260) state in B 
meson decay, mediated by an effective vertex operator Oi. 
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The process B —> Y (4260 )K occurs via weak decay of 
the b quark, while the u quark is a spectator. The Y 
meson as a mixed state of tetraquark and charmonium 
interacts via cc component of the weak current. In effec¬ 
tive theory, at the scale p ~ mb <C mw, the weak decay 
is treated as a four-quark local interaction described by 
the effective Hamiltonian (see Fig. ©: 


% = ~^ VcbV * t 


c 2 {ii) + ) o 2 


, ( 2 ) 


meson pole Q 2 = -Wy, u sing the QCDSR approach for 
the three-point correlator [301 ]: 

n„(p,p') =y d 4 xd 4 ye^'- x - p -y\0\T{J^(0) x 

x J K (x)Jl(y) }|0), (9) 

where the weak current, , is defined in (0 and the in¬ 
terpolating currents of the B and K pseudoscalar mesons 
are: 


where Vik are CKM matrix elements, Ci(/x) and C 2 (y) 
are short distance Wilson coefficients computed at the 
renormalization scale p, ~ 0(mb). The four-quark effec¬ 
tive operator is 0 2 = , with 

jJ = aC„b, jf)=cT M c, (3) 

and T^ = 7 m (1 - 75 ). 

Using factorization, the decay amplitude of the process 
is calculated from the Hamiltonian m, by splitting the 
matrix element in two pieces: 

m = ‘^vm(c, + f) 

x (B(p)\jW\K(p'))(Y( q )\J^\0), (4) 

where p = p'+q. Following Ref. [29l ]. the matrix elements 
in Eq. © are parametrized as: 


Jk — ^ ^a75 1 JB — f • (10) 

The obtained result for the form factor was f2j|: 


2 = (17.55 ±0.04) GeV 2 

/+lW ' (105.0 ± 1.8) GeV 2 + Q 2 ' 


( 11 ) 


For the decay width calculation, we need the value of 
the form factor at Q 2 = —my, where my is the mass of 
the V(4260) meson. Using my = (4251 ±9) MeV [jd! we 
get: 

/+(Q 2 )| Q 2 = _ m 2 , = 0.206 ± 0.004. ( 12 ) 

The parameter Xw can also be determined using the 
QCDSR approach for the two-point correlator: 


Iw«) =ijd 4 y e^(0|T{jJ(y)4 5 c) (0)}|0), (13) 


(Y(q)\J^\0) = X w e;(q), 

and 


(5) 


where the current Jis defined in ([3]). For the Y meson 
we will follow [ 2 (il] and consider a mixed charmonium- 
tetraquark current: 


(B(p)\jW\K(p')) = f + {q 2 )(p, + */„) + f-(q 2 )(p,-p'J . 

( 6 ) 

The parameter A w in © gives the coupling between the 
current J j l cc ' > and the Y state. The form factors f±(q 2 ) 
describe the weak transition B — ► K. Hence we can see 
that the factorization of the matrix element describes the 
decay as two separated sub-processes. 

The decay width for the process B~ —> Y (4260)Rf - is 
given by 

T(B -j. YK) = T ^ r y/x(m 2 B ,m 2 K ,m 2 Y ), (7) 

0 B 

with X(x, y , z ) = x 2 + y 2 + z 2 — 2 xy — 2 xz — 2 yz. The 
invariant amplitude squared can be obtained from ©, 
using © and ©: 

= w 2 4 + t) 2 


= sin 9 J(f 1 + cos 6 , 


where 


^abc^dec 


/(4) = 

M V2 


{Qa C 'y5Cb)(q d 'r^5Cc^) 


(q a C'757/xCh)((?d75C'c e ) 


4 2) = ^( 99 ) (ca7/.Ca) = J } 2) 


V 2 ' 


(14) 


(15) 

(16) 


In Eq. (EM is the mixing angle that was determined in 
[H to be: 6 = (53.0 ±0.5)°. 

Inserting the currents © and (TTH) in the correlator we 
have in the OPE side of the sum rule 

n °„ E (q) = sin9H^ 2 (q) ± ^ cos6»n 2 ’ 2 (g), (17) 

where 


x X{m 2 B ,m 2 K ,mlr)X 2 w fl. 


( 8 ) 


The coupling constant /+ was determined in Ref. [291] 
through extrapolation of the form factor f + (Q 2 ) to the 


n^(<?) = i J d A y e^{0\T{jW(y)J vlSc) {0)}\0) 
n 2,2 (?) =ijd 4 y e iqy (0\T{j'( 2 \y)J v{ , c) (0)}\0) .(18) 
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TABLE I. QCD input parameters. 


Parameters 

Values 

m c 

(1.23 - 1.47) GeV 

(w) 

— (0.23 ± 0.03) 3 GeV 3 

(dsG 2 ) 

(0.88 ±0.25) GeV 4 

ml = (qGq)/(qq) 

(0.8 ±0.1) GeV 2 


(cc) 

Only the vector part of the current Jc, ' contributes to 
the correlators in Eq. m • Therefore, these correlators 
are the same as the ones calculated in Ref. [26| for the 
mass of the Y (4260). 

To evaluate the phenomenological side we insert inter¬ 
mediate states of the Y: 


<T(<z) = 


9 2 

q z — rriy 


(0\J^\Y(q))(Y(q)\J^\0) , 


Q 2 


9 




°Y 


(19) 


where q 2 = —Q 2 , and we have used the definition 0 and 


(0\J^\Y(q))=X Y € li (q). ( 20 ) 


The parameter Ay, that defines the coupling between the 
current J 4 and the Y meson, was determined in Ref. j26] 
to be: Ay = (2.00 ± 0.23) x 10" 2 GeV 5 . 

As usual in the QCDSR approach, we perform a Borel 
transform to Q 2 —> M B to improve the matching between 
both sides of the sum rules. After performing the Borel 
transform in both sides of the sum rule we get in the g 
structure: 


Aw Aye = ^n 4 > 2 (M|) + -^-cos6»n 2 ’ 2 (M|)(21) 
V2 \/2 


where the invariant functions II 2 ’ 2 (Mg) and II 4 ’ 2 (M|>) 
are written in terms of a dispersion relation, 


OO 

n(M 2 ) = J ds e -s / M B p( s ) , (22) 

4777.2 


with their respective spectral densities p 2,2 (s) and p 4,2 (s) 
given in Appendix. 

We perform the calculation of the coupling parame¬ 
ter A w using the same values for the masses and QCD 
condensates as in Ref. [26] which are listed in Table 
H] To be consistent with the calculation of Ay we also 
use the same region in the threshold parameter sq as in 
Ref. (26|: = (4-70 ± 0.10) GeV. As one can see in 

Fig. 1 the region where we get M^-stability is given by: 
(8.0 <M 2 b < 25.0) GeV 2 . 

Taking into account the variation in the Borel mass 
parameter, in the continuum threshold, in the quark con¬ 
densate, in the coupling constant Ay and in the mixing 
angle 6, the result for the A w parameter is: 

Aw = (0.90 ± 0.32) GeV 2 . 


Thus we can calculate the decay width in Eq. 0 by 
using the values of /+(— M 2 ) and A w, determined in 
Eqs. m and (l23l) . The branching ratio is evaluated 
dividing the result by the total width of the B meson 
T to t = 4.280 x 10 " 4 eV: 

B(B -Y Y(4260)iT) = (1.34 ± 0.47) x 10 “ 6 , (24) 

where we have used the CKM parameters 14s = 1-023, 
14b = 40.6 x 10 -3 ,3l|, and the Wilson coefficients 
Ci{p) = 1.082, C 2 (/x) = —0.185, computed at p, = mb 
and A M s = 225 MeV H3- 

In order to compare the branching ratio in Eq. (l24l) 
with the branching fraction obtained experimentally in 
Eq. 0, we might use the results found in Ref. (26j : 

Z?(Y(4260) —$■ J/ip 7 T+ 7 T-) = (4.3 ± 0.9) x 10 " 2 , (25) 

and then, considering the uncertainties, we can estimate 
B Y > 3.0 x 10~ 8 . However, it is important to notice 
that the authors in Ref. [26[ have considered two pions 
in the final state coming only from intermediate states, 
e.g. a and /o(980) mesons, which could indicate that the 
result in Eq. (E5l) can be underestimated. In this sense, 
considering that the main decay channel observed for the 
Y(4260) state is into J/ij} 7 r + 7 r _ , we would naively expect 
that the branching ratio into this channel could also be 
£>(Y(4260) — > J/^ 7 r + 7 r - ) ~ 1.0, which would lead to the 
following result, B Y < 1.8 x 10^ 6 . Therefore, we obtain 
an interval on the branching fraction 

3.0 x 10 " 8 < B Y < 1.8 x 10 ~ 6 (26) 

which is in agreement with the experimental upper limit 
reported by Babar Collaboration given in Eq. 0. In 
general the experimental evaluation of the branching 
fraction takes into account additional factors related 
to the numbers of reconstructed events for the final 
state (J/i/j 7 r + 7 r - K ), for the reference process (B — > 
Y(4260) K ), and for the respective reconstruction effi¬ 
ciencies. However, since such information has not been 



FIG. 2. The coupling parameter Aw as a function of Mg, for 
different values of the continuum threshold. 


(23) 
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provided in Ref. [|], we have neglected these factors in 
the calculation of the branching fraction By. Therefore, 
the comparison of our result with the experimental result 
could be affected by these differences. 

In conclusion, we have used the QCDSR approach to 
evaluate the production of the Y (4260) state, considered 
as a mixed charmonium-tetraquark state, in the decay 
B —y YK. Using the factorization hypothesis, we find 
that the sum rules result in Eq. (1241) . is compatible with 
the experimental upper limit. Our result can be inter¬ 
preted as a lower limit for the branching ratio, since we 
did not considered the non-factorizable contributions. 

Our result was obtained by considering the mixing an¬ 
gle in Eq. (1141) in the range 0 = (53.0 ± 0.5)°. This angle 
was determined in Ref. [26| where the mass and the decay 
width of the F(4260) in the channel were de¬ 

termined in agreement with experimental values. There¬ 
fore, since there is no new free parameter in the present 
analysis, the result presented here strengthens the con¬ 
clusion reached in |26| that the T(4260) is probably a 
mixture between a cc state and a tetraquark state. 

As discussed in (29j, it is not simple to determine the 
charmonium and the tetraquark contribution to the state 
described by the current in Eq. m- From Eq. (ED one 
can see that, besides the sin 0 , the cc component of the 
current is multiplied by a dimensional parameter, the 
quark condensate, in order to have the same dimension 
of the tetraquark part of the current. Therefore, it is 
not clear that only the angle in Eq. ED determines the 
percentage of each component. One possible way to eval¬ 
uate the importance of each part of the current it is to 
analyze what one would get for the production rate with 
each component, i.e., using 9 = 0 and 90° in Eq. ED- 
Doing this we get respectively for the pure tetraquark 
and pure charmonium: 

B{B -4 Y tetra I<) = (1.25 ± 0.23) x 10 " 6 , (27) 

B{B Y 5c K) = (1.14 ± 0.20) x 10 “ 5 . (28) 

Comparing the results for the pure states with the one 
for the mixed state (EDi we can see that the branching 
ratio for the pure tetraquark is one order smaller, while 
the pure charmonium is larger. From these results we see 
that the cc part of the state plays a very important role 
in the determination of the branching ratio. On the other 
hand, in the decay Y J/^ 7 T + 7 r _ , the width obtained 
in our approach for a pure cc state is [26j : 

r(Yec -4 J/V>tttt) = 0 , (29) 


the T(4260), we may say that both components are ex¬ 
tremely important, and that, in our approach, it is not 
possible to explain all the experimental data about the 
T(4260) with only one component. 
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Appendix A: Spectral Densities for the Two-point 
Correlation Function 


We list the spectral densities for the invariant functions 
related to the coupling between the current and the 
Y (4260) state. We consider the OPE contributions up to 
dimension-five condensates and keep terms at leading or¬ 
der in a s . In order to retain the heavy quark mass finite, 
we use the momentum-space expression for the heavy 
quark propagator. We calculate the light quark part 
of the correlation function in the coordinate-space and 
use the Schwinger parametrization to evaluate the heavy 
quark part of the correlator. For the d 4 y integration in 
Eq. (1131) , we use again the Schwinger parametrization, af¬ 
ter a Wick rotation. Finally, the result of these integrals 
are given in terms of logarithmic functions through which 
we extract the spectral densities. The same technique can 
be used for evaluating the condensate contributions. 

Then, in the g structure, we evaluate the spectral 
densities for the II 2 ’ 2 (M B ) function, 


2,2 ^ ^ _ rrii A, , l\ , ( g 2 s G 2 ) v l\ 

P {s) -4^ V \ 2+ x) + ^ 8^ Ml 


M 2 b x 


(“-§) 


\M%x 


B L 

2 'A—i 

X 


)(> 


,(A1) 


and for the II 2 ’ 4 (M|) function, 


P IA {s) 


ml (gg) 

12t r 2 



(qGg) 

247T 2 




and, therefore, the tetraquark part of the state is the only 
one that contributes to this decay, playing an essential 
role in the determination of this decay width. 

Therefore, although we can not determine the per¬ 
centages of the cc and the tetraquark components in 


where we have used the definitions 

x = m 2 c /s (A3) 

v = Vl — 4a; . (A4) 


[1] B. Aubert et al. [BaBar Collaboration], Phys. Rev. Lett. [2] Q. He et al. [CLEO Collaboration], Phys. Rev. D 74, 
95, 142001 (2005). 












5 


091104(R) (2006); C.Z. Yuan et al. [Belle Collaboration], 
Phys. Rev. Lett. 99, 182004 (2007). 

[3] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 73, 
011101 (2006). 

[4] S. L. Zhu, Int. J. Mod. Phys. E 17, 283 (2008) 
hep-ph/0703225 . 

[5] M. Nielsen, F. S. Navarra and S. H. Lee, Phys. Rept. 
497, 41 (2010) _arXiv:0911.1958 . 

[6] N. Brambilla, et al., Eur. Phys. J. C71, 1534 (2011) 
arXiv: 1010.5827. 

[7] M. Nielsen and F. S. Navarra, Mod. Phys. Lett. A 29, 
1430005 (2014) arXiv:1401.2913 . 

[8] M. Ablikim et al. [BESIII Collaboration], Phys. Rev. 
Lett. 110 , 252001 (2013). 

[9] Z.Q. Liu et al. [BELLE Collaboration], Phys. Rev. Lett. 
110, 252002 (2013). 

[10] T. Xiao, S. Dobbs, A. Tomaradze and K.K. Seth, Phys. 
Lett. B 727, 366 (2013). 

[11] L. Maiani, V. Riquer, F. Piccinini and A. D. Polosa, 
Phys. Rev. D72, 031502 (2005). 

[12] G. J. Ding, Phys. Rev. D79, 014001 (2009). 

[13] C. Z. Yuan, P. Wang and X. H. Mo, Phys. Lett. B634, 
399 (2006). 

[14] X. Liu, X. Q. Zeng and X. Q. Li, Phys. Rev. D72, 054023 
(2005). 

[15] A. Martinez Torres, K. P. Khemchandani, D. Gamer- 
mann, E. Oset, Phys. Rev. D80, 094012 (2009) 
arXiv:0906.5333 . 

[16] S. L. Zhu, Phys. Lett. B625, 212 (2005). 

[17] C. F. Qiao, Phys. Lett. B639, 263 (2006). 

[18] E. van Beveren and G. Rupp, arXiv:hep-ph/0605317 

[19] E. van Beveren and G. Rupp, arXiv:0904.4351 

[20] E. van Beveren and G. Rupp, Phys. Rev. D79, 111501 


(2009). 

[21] M.A. Shifman, A.I. and Vainshtein and V.I. Zakharov, 
Nucl. Phys. B 147, 385 (1979). 

[22] L.J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 
127, 1 (1985). 

[23] For a review and references to original works, see e.g., 
S. Narison, QCD as a theory of hadrons, Cambridge 
Monogr. Part. Phys. Nucl. Phys. Cosmol. 17, 1 (2002) 
hep-h/0205006 ; QCD spectral sum rules , World Sci. 
Led. Notes Phys. 26, 1 (1989); Acta Phys. Pol. B 26, 
687 (1995); Riv. Nuov. Cim. 10N2, 1 (1987); Phys. 
Rept. 84, 263 (1982). 

[24] R.M. Albuquerque and M. Nielsen, Nucl. Phys. A815, 
53 (2009); Erratum-ibid. A857 (2011) 48. 

[25] R. M. Albuquerque, M. Nielsen and R. R. da Silva, Phys. 
Rev. D 84, 116004 (2011) arXiv:1110.2113 [hep-ph]]. 

[26] J. M. Dias, R. M. Albuquerque, M. Nielsen and 

C. M. Zanetti, Phys. Rev. D 86, 116012 (2012) 

arXiv: 1209.6592 . 

[27] R. D. Matheus, F. S. Navarra, M. Nielsen and 
C. M. Zanetti, Phys. Rev. D 80, 056002 (2009) 
arXiv:0907.2683 . 

[28] M. Nielsen and C. M. Zanetti, Phys. Rev. D 82, 116002 
(2010) arXiv: 1006.0467 [hep-ph]]. 

[29] C. M. Zanetti, M. Nielsen and R. D. Matheus, Phys. Lett. 
B 702, 359 (2011) arXiv: 1105.1343 [hep-ph]]. 

[30] M. E. Bracco, M. Chiapparini, F. S. Navarra and 
M. Nielsen, Prog. Part. Nucl. Phys. 67, 1019 (2012) 
arXiv:1104.2864 . 

[31] K.A. Olive et al. (Particle Data Group), Chin. Phys. C 
38, 090001 (2014). 

[32] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. 
Mod. Phys. 68, 1125 (1996) arXiv:hep-ph/9512380 . 




